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Dust in the early Universe: Evidence for non-stellar dust production 
or observational errors? 
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ABSTRACT 

Observations have revealed unexpectedly large amounts of dust in high-redshift galaxies and 
its origin is still much debated. Valiante et al. (2009, MNRAS, 397, 1661) suggested the 
net stellar dust production of the quasar host galaxy SDSS Jl 148+5251 may be sufficient to 
explain the large dust mass detected in this galaxy, albeit under some very special assumptions 
(e.g., 'closed box' evolution and a rather high gas mass). Here it is shown that since accretion 
of essentially pristine material may lower the efficiency of dust formation significantly, and 
the observationally derived dust-to-gas ratios for these high-redshift galaxies are remarkably 
high, stellar dust production is likely insufficient. A model including metallicity-dependent, 
non-stellar dust formation ('secondary dust' ) is presented. The required contribution from this 
non-stellar dust component appears too large, however. If all observational constraints are 
to be met, the resultant dust-to-metals ratio is close to unity, which means that almost all 
interstellar metals exist in the form dust. This is a very unlikely situation and suggests the 
large dust-to-gas ratios at high-redshifts may be due to observational uncertainties and/or or 
incorrect calibration of conversion factors for gas and dust tracers. 

Key words: Galaxies: evolution, high-redshift, starburst; Quasars: general, individual: SDSS 
Jl 148+5251; Stars: AGB and post-AGB, supernovae: general; ISM: dust 



1 INTRODUCTION 

Observations of damped Lya systems (DLAs), quasars and gamma- 
ray bursts (GRBs) have shown that relatively l arge amounts 
of dust are present in high-z galaxies (see , e.g. [L edoux et al. 
' | 2002l: iBertoldi et al.1120031; iRobson et alj|2004 iBeelen et al.ll200rl 

iMichalowski et alj|200 8. 2010a b). This raises questions about the 
origins of cosmic dust and its effect on cosmological observations. 
Most of the dust present in the Milky Way today is probably pro- 
duced in the envelopes of evolved (age > 1 Gyr), low-mass stars. 
Observations of large amounts of dust present in high-z galaxies 
with ages of less than 1 Gyr, on the other hand, suggest that even 
though low-mass stars seem to be the dominant dust producers 
in the present universe, this might not alw ays have been the case 
jBertoldi et alj2003l : lMaiolino et alj2004h . 

Core-collapse supernovae (SNe) arising from early genera- 
tions of stars have been sugge sted as the main sources of dust 
in the early Universe dTodini & Ferrarall200ll: iNozawa et al .1 120031 ; 

iDwek et alj|2007l : iBianchi & Schneiderll2OO70 Models of dust for- 
matio n in core-collapse SNe dTodini & Ferrara 2 00lllNozawa et all 
2008) predict large enough amounts of dust to account for the dust 
seen at high z. Direct observational evidence for core-collapse SNe 
as a major dust pr oducers is however still inconclusive, even for 
the local Universe dMeikle et alj2007l) . Observational results seem 
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to suggest the onset of dust condensation in SN ejecta, although 
the SN contribution to total mass of cosmic dust appears to be mi - 
nor in most cases (see, e.g.. iKotak et~afll200d [20091 : lGallll2010h . 
However, a few intersting cases where significant am ounts of dust 
are detected do exist, su ch as Keplers supe rnova I Gomez et al.l 
20091: iMorgan et alj|2003h and Cassiopeia A dBarlow et alj|2010l : 



Dunne et"al]| 20091) . There is also evi dence for a significa nt amount 
dust (0.4M o ) surrounding Eta Carina dGomez et al.l2010h . Whether 
the dust present in the Kepler and Cas A SN remnants is produced 
by the SNe themselves or the dust is formed from metals in the 
surrounding ISM is unclear. It is also uncertain if the dust which 
can be formed in the winds of massive stars (such as Eta Carina) 
survive the shock wave of the subsequent SN explosion. However, 
it is not possible to rule out SNe as significant dust producers, al- 
though a recent analysis of the existing observational constraints 
on dust productivity in SNe has shown that, on average, SNe seem 
to be converting a much smaller fraction of their metals into dust 
(correspon ding to a fe w times 10~ 4 -10~ 2 M o of dust) compared to 
AGB stars dGalfcOlOh . 

Assum ing SNe as the only source of dust in the early Universe, 
IDwek et all ( 2007) suggest a dust yield o f 1M G of dust per SN to ac- 
count for dust masses in high-z quasars. IValiante et al.ld2009l) used 
the SN dust yields including dust destruct ion by the revese shock 
computed by IBianchi & Schneidei] d2007l) and the currently only 
existing sets of dust yield s for lo w an d intermediate mass (LIM) 
stars by iFerrarotti & Gaill (§006) and [Zhukovska, Gail & Trieioffl 
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J2008I) in combination with the star-formation history obtained by 
simulati on of a hierarch ical assembly in a canonical ACDM cos- 
mology jLietal]l2007h for the host galaxy of the second most 
distant (z = 6.42) quasar ever detected: SDSS J1148+5251. This 
galaxy appears to contain an unexpectedly high dust mass (given 
its redshift) which has raised doubts about w hether stars alone can 
produce so much dust in such a short time. IValiante et al.l J2009I) 
argued that stellar sources can account for the dust detected by ob- 
servations and show also that the contribution of AGB stars to dust 
production in the early Universe m ay be far from negligible. 

R ecently. IPipino et alj 201 lh reached a conclusion similar to 
that of IValiante et al.l d2009 ). although their model includes also 
nonste llar sources of dus t. However, t he mo del by IValiante et al.l 
(2009) as well as that of IPipino et al.l d201ll) assumes a gas mass 
that is much larger (about an order of magnitude) than the observed 
mass of molecular gas of ~ 10 10 M G (Walt er et al.ll2004h . A more 
gas rich galax y will obvio usly be able to hold a greater mass of 
dust, which led Gall (2010) to suggest that the gas content of many 
dust rich galaxies at high redshifts may be greatly underestimated 
(their preferred models hasafew times 10 U M of gas after 1 Gyr). 

An intriguing discovery which has complicated the picture is 
that not all quasar host galaxies at high redshifts show significant 
dust compon ents. In fact, som e of them may be essentially dust 
free (see, e.g. Jjiang et al.l200or) . Quasars at high redshifts are pow- 
ered by supermas s ive black holes w ith masses > 10 9 M o (see, e.g., 
IVeste rgaard 2o61 |jiang etafl l2006). which may be part of the ex- 
planatio n to why some qua sar hosts are very dust rich while others 
are not. Ijiang et all d201 Ol) have shown that there is a correlation 
between the estimated size of the central black hole and the derived 
dust masses. This finding may lend some support to a hypothesis 
presented by lElvis et al. I d2002h . that dust may form in significant 
quantities in the outflows of quasars. 

In this paper it will be shown that the stellar dust production 
required to explain the observed dust mass in SDSS Jl 148+5251 
and in other high-z galaxies is higher than what appears reasonable 
taking both theoretical and observational constraints into account. 
Is some kind of 'secondary' dust production, possibly associated 
with a quasar outflow, required to explain the detections of large 
dust masses? Or are the high dust-to-gas ratios at high redshifts 
simply an artefact due to incorrect calibration of conversion factors 
for gas and dust tracers? 



2 THEORY AND MODEL 

2.1 Cosmology and time of galaxy formation 

Adopting a canonical ACDM cosmology, a closed form express ion 
for the expansion scale factor 5 (f) can be obtained (|Gr 0nll2OO2l) , 
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where Hq is Hubble's constant, Q. M is the fractional matter energy 
density today, and Q A is the fractional dark energy density today. 
Combining this expression with z + 1 = S(to)/S(t) gives the time- 
redshift relation 

ft 

(z+ir 3/2 



7-1 



2 H- 
W) = T 77777 smh 



(2) 



Here, a standard scenario with H = 70 km s _I Mpc~\ fl M = 0.3 
and Q A = 0.7 is assumed. Galaxy formation is assumed to start 
at z = 14, which corresponds to f a » 290 Myr. This choice of 



redshift is somewhat arbitrary, but the first galaxies are believed to 
have formed during the reionisation epoch, i.e., between z = 20 and 
z = 6. Galaxy formation at z = 14 may thus be representative for 
when the vast majority of galaxies formed in the early Universe. 

According to the current cosmological paradigm (see Sect. 
12. IK galaxies are formed by mergers with other galaxies and/or 
infall of baryons at a certain rate M\ . The time scale of infall (mass 
assembly) may be defined as r, = M g /M[. It is reasonable to as- 
sume the star-formation rate somehow reflects the rate of infall. In 
particular, it is especially convenient (mathematically) to assume 
M{ = aM*, which is equival ent to assumi ng that the gas mass is 
constant with respect to time (Larson 1972). This simplistic model 
of the rate of infall will be used throughout this paper. 



2.2 Dust destruction 

Dust in the ISM is expected to be destroyed by kinetic-energy injec- 
tion to the ISM by SNe, but dust grains may also grow in the ISM, 
in particular, the outer parts of SN remnants. Further dust growth 
may also take place in molecular clouds. However, dust destruction 
in the ISM is probably dominating over nucleation (which rarely 
happens in the ISM), and it is thus fair to assume that introducing 
an effective rate of dust destruction in the ISM Mism i s suffic ient 
in order to model the dust cycle. Following Dwe k et all d2007h the 
dust destruction time-scale is 



misM ^s 



(3) 



where M g is the gas mass, m ISM is the effective gas mass cleared 
of dust by each SN event, and i?sN is the SN rate, which may be 
approximated as 



#snC) * M 



*(r) I <p(m)dm 
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where m n is the mass of the most massive SN-progenitors. The in- 
tegral in Eq. (|4) is a constant with repspect to time, hence the time 
scale Td may be expressed as 



M g 



(5) 



where (5ism will be referred to as the efficiency of dust destruction in 
the ISM for each generation of stars. For a normal IMF and misM ~ 
IOOOMq (the preferred value for the Milky Way, see Dw ek et all 
120071 and references therein), this parameter is (5ism ~ 10.0. 

When considering the contribution from AGB stars there is 
a time-lag due to the lifetimes of these stars and therfore the de- 
struction of dust depends on how the SN rate evolves. If the rate of 
star formation has decreased sufficienly when the release of new 
dust from AGB stars happens, the dust destruction by SNe will 
be much less due to the lower SN rate. Effectively the parame- 
ter S\sm may have to be smaller than the number estimated above. 
Dwe k et al T j2007l) considered mi S M a free parameter, and found 
that with misM = 100M o and a top-heavy IMF (shallower than the 
ISalpeterl 19551 IMF) the estimated dust mass of SDSS Jl 148+525 1 
could be explained in terms of dust production in SNe/high-mass 
stars. If misM = 100M o is adopted as the effective gas mass cleared 
of dust by a SN-event, then 6\sm ~ L0. 

Depletion of dust in the ISM need not be due to any actual 
destruction of dust. If a galactic wind is present it may have es- 
sentially the same effect, although gas and atomic metals will be 
lost too. Since a galactic wind is beleived to be mainly driven by 
the kinetic-energy injection by SNe, it is fair to assume that M w is 
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Figure 1. Left: theoretical remnant masses as function of stellar mass compared to a simple power-law. Right: lock-up fraction as function of the turn-over 
mass m c and the lower mass limit m T for stars contributing to the build-up of dust and metals. 



proportional to the star-formation rate M t (see Eq. |4j. Hence, the 
outflow time-scale of the wind is 



(6) 



where, by analogy with Eq. l|5}, M w is the galactic-wind mass loss, 
and <5 W is the galactic 'wind efficiency'. However, for simplicity, 
only actual destruction of dust will be considered in the following, 
but the reader should bear in mind that a galactic wind may have a 
similar effect on the dust-to-gas ratio Zj. 



where M t is the star-formation rate a is the stellar lock-up fraction, 
i.e., the fraction (typically about 60-80%) of the gas mass converted 
into stars that will remain locked up in stellar remnants and very 
long-lived stars, y$ is the dust yield for a single generation of stars 
and Mi is the infall rate. Combining the equations above and the 
prescriptions for infall and dust destruction one arrives at 



dZ A dM s l6 lSM a \dM 
Ms ~ =>' J — "(— Zd 



(10) 



where Zd = M A /M„ is the dust-to-gas mass ratio and {3 = 0, 1 for 
the closed-box and extreme infall cases, respectively. 



2.3 Equations of dust evolution 

In order to obtain analytical solutions of adequate simplicity, the 
instantaneous recycling approximation (IRA, which esse ntially 
mean s all stars are assumed to have negligible lifetimes, see lPagej 
1 19971) is employed throughout the rest of this paper. No delayed el- 
ement production due to stellar life-times is considered. This may, 
or may not, be a good assumption, but since the purpose of this 
paper is to provide 'proof of concept' and to study the dust compo- 
nent primarely during earliest phase of galaxy evolution, the IRA 
is a justified simplification. 

A high-z galaxy is assumed to be adequately modelled by a 
one-zone model and the production of metals have a negligible ef- 
fect on the evolution of the gas mass. Including infall of pristine 
gas and dust destruction in the ISM the eqautions for the evolution 
of stars, gas and dust then becomes 



dM s 
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dM g 
~dF 
dM A 
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2.4 Lock-up fraction and stellar dust yields 

The primar)[] dust yield Vd for a single generation of stars is defined 



i r*» 

a J,„ r 



Pd(m) m (f)(m) dm. 



(11) 



In Eq. dl lb above, p A is the fraction of the initial mass m of a star 
ejected in the form of newly produced dust and <f>(m) is the mass- 
normalised IMF. A convenient functional form of the IMF which 
in cludes the turn -down at low stellar masses is the form suggested 
bv lLarsod ( fl998l) . 



(p(m) = 0o m (1+Jr) exp ( ) , 

\ ml 



(12) 



where m c is the characteristic mass around which the turn-down be- 
gins and (p is a normalisation constant obtained from the condition 



m (j>(m) dm = 1 . 



(13) 



The upper and lower mass-cuts, m u and tn\, respectively, represents 
the high-mass trunction and the brown-dwarf limit. 

1 Dust is here assumed to form from newly synthesised elements. 
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Table 1. IMF parameters and the resultant stellar lock-up fractions for the 
two considered IMFs. 
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[M Q ] 


[M Q ] 


[M G ] 


[M Q ] 




Normal 


1.35 


0.35 


0.1 


100.0 


3.0 


0.63 


Top-heavy 


1.35 


10.0 


0.1 


100.0 


3.0 


0.10 



In the IRA the gas-return fraction R of a generation of stars is 
obtained from the convolution of the IMF and the stellar remnant 
mass w(m), i.e., 



*= r 



[m — w(m)] <p(m) dm. 



(14) 



The remnant mass w(m) is a function of the initial stellar mass 
m, which may be approximated by a simple power-law, w(m) = 
k m", where k is a constant. Usin g the remnant masses obtained 
bv Ivan den H oek & Groenewegen] dl997h and IWooslev & Weaverl 
dl995h for LIM and high mass (HM) stars, respectively, one ob- 
tains k = 0.5 and n = 0.4 for stellar masses less than m = 35M Q 
(see Fig.[T]l. 

The lock-up fraction a, used in the model equations discussed 
in Sect. 12. 31 is obtained as a = 1 - R. Adopting an IMF of the form 
shown above in Eq. J22b and defining G(a, b, c, z) = T(a, z/c) - 
T(a, z/b), where T(a, z) is the incomplete gamma-function, one ar- 
rives at the expression (which holds for x > 1, n < 1) 

G(x' ,m a ,m T ,m c ) - km" G(x' - n' ,m a ,m T ,m c ) 

R = 7T~, \ ' (I 5 ) 

G(x' ,m u ,m\,m c ) 

where xf = x — 1, n' = n - 1 and m T is the minimum stellar mass 
effectively contributing to the matter cycle at the considered red- 
shift. With the parameter values given in Table Q~]and m T = 3.0M o 
(the initial mass for which the stellar lifetime corresponds to the 
typi cal age o f gala xies at z = 5 - 6) the lock-up fraction for a nor- 
mal |Lareon| d 19981) IMF (m c = 0.35M G ) becomes a = 0.63, and 
in the top-heavy (m c = 10.0M o ) case a = 0.10. As shown in the 
right panel of Fig.Q] a is a steep function of the turn-over mass m c , 
but only weakly dependent on m T . Hence, the choice of m T is not 
critic al fo r either a, nor the effective primary yield y d . 

computed the total dust productivity of a genera- 
tion of stars considering different IMFs and st ellar dust-pr oduction 
efficiencies, e d (m) = y d (m)/y z (m). Accordig to lGalll J2010h the dust 
production efficiency derived from observations of SN remnants 
is well-approximated by a power law, which con nects nicely with 
theor etically expected efficiencies for AGB stars l iFerrarotti & Gaill 
2006). If this is a general law for stellar dust production, SNe 
are insignificant dust producers when weighted by an IMF. The- 
oretical work, as well as a handfull of observational studies, al- 
lows for a much highe r efficiency in SNe. Using the results of 
iTodini & Ferraral J200 ll) . iGalll d20 1 0b find that AGB stars would be 
insignificant relative to SNe as dust producers if no dust destruction 
of SN-formed dust occur due to the reverse shock. 

The integrated (IMF-weighted) yields calculated by IGalll 

(2010) is used in this paper, i.e., three cases are consid- 
ered (see Table |2j: 'observed' SN dust yields, theoretical SN 
yields with dust destructio n due to reverse shocks ('high' yield, 
iBianchi & Schneider||2007[). and t he theoretical 'upper limit' ac- 
cording to Todi ni & Ferraral teooil referred to as 'maximal' yield) 
- all three in combination with the theoretical results for A GB stars 
dFerrarotti & Gaill2006t Izhukovska, Gail & Trieloflll2008l) . 



Table 2. Yields used in this paper. 



IMF 



jnax 
" d 



Normal 
Top-heavy 



8.81 ■ 10~ 4 
3.51 ■ 10- 3 



i.52 • io-3 i.08 -nr 2 

1.55 -10- 2 1.88 Kr 1 



2.5 Evolution of the dust-to-gas ratio 

2.5.1 Closed-box evolution 

A 'closed-box' evolution is obviously incompatible with the sce- 
nario described above in Sect. 12.1 1 But fo r completenes, and si nce 
it corre sponds to the scen ario studied by IValiante et al.l d2009h as 
well as iGall et alj J2010b1) . it is included also here. If there is no 
infall, the rate of change of stellar mass is 



(16) 



1 = aM* = 

dt dt 

Assuming no dust destruction, combination of Eq. d!6t and Eq. J 10b 
gives the classical closed-box solution, 



• y d In 1 1 + 



M„ 



(17) 



Replacing y d with the metal yield y z , this is also the solution for the 
metallicity Z in a closed-box s cenario. If the dust-destruction term 
is included, the solution reads (Edmunds 20 0ll : lDwek et al.ll2007l) . 



v 



M s 

1 - 1 + — 



dlSM 

a 



(18) 



This solution differs from the classical closed-box soultion in that 
■Zd — > yi/v as M s /M g — » oo. Assuming that dust grains may form 
and grow out of metals produced by previous generations of stars, 
i.e., a kind of 'secondary' dust, one may redefine the dust yield as 



M0 = yd + eZ(t), 



(19) 



where e is the fraction of the metals present in the ISM that will 
form dust grains. There is a natural upper limit to the value of e 
since Z d < Z, i.e., the dust mass cannot exceed the total mass of 
metals available. 

With the modifi cation above, the solution for Z d becomes (see 
also lEdmundsll200ll) 



Vd 



yz 



Z, = — 1 1 - w — 

yd 



i-i + 



M s 



+ u>y 7 In 1 H 

1 AfJ 



(20) 



where w = e/v. For the special case Sism = (no dust destruction) 
the solution is 

2 

(21) 



MA e 
Z d=yd ln(l + - + -y z 



M„ 



2.5.2 Galaxy formation by infall 

Assume the rate of gas consumption is exactly balanced by the rate 
of infall. With this assumption made, and no dust destructio n term, 
the so lution to Eq. d 1 0b is the 'extreme infall model' by iLarsonl 
dl972h . i.e., 



yd 



exp 



M s 

"mT 



(22) 



Adding dust destru ction, the solu tion for Z d has the same mathe- 
matical form as the |Earson| Jl972|) solution, but with an extra pa- 
rameter y, 
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7 yi 



exp|- 7 — 



7 = v + 1, 



(23) 



which reduces to Eq. d22t if (5 1SM = (7 = 1). Using the modified 
yield introduced in Eq. d!9t the solution for Zd becomes 



— I 

y 



1 — OJ 



1 - exp -7 



'~mI 



+aiy z 



1 - exp 



,(24) 



where at is defined as before, for the closed-box case. For the spe- 
cial case (Sism = (no dust destruction) the solution is 



1 Ms 



+ eyz 



1 ~~ ex P ~ 77" 1 + 77" 



M„ 



Mo 



(25) 



2.5.3 Evolved systems 

Even high-z galaxies may have quite evolved stellar populations, 
e.g., they may be in (or approaching) a 'post-starburst state' after 
an intense initial episode of star formation. One should note that 
the dust-to-gas ratio Zd in the solutions above can never decrease 
(although Md can) since dM s /dt is always positive. Hence, it is of 
interest to know whether Z d is approaching a constant value and 
thus has an upper limit or not. 

For a closed-box scenario without dust destruction by SNe 
(<5ism = 0), Z d — > 00 as M s /M„ — » 00. Including dust destruc- 
tion (<5 ISM ± 0), there is an upper limit since then Z d — > Vd/v as 
M s /M g — > 00, given that y& is treated as constant. Including 'sec- 
ondary' dust production, Z d — » 00 since the second terms of Eqs. 
j20\ and d21t are simple functions the closed-box solution without 
dust destruction, which are monotonically increasing. 

Any infall model where the rate of infa l l is pr oportional to 
the star-formation rate - not only the iLarsonl (Il972b model - has 
always an upper limit to Zd. In Larson's extreme infall model the 
metallicity can never exceed the yield. A similar upper limit exists 
also for a dust evolution model with infall, although the limit may 
be lower due to dust destruction. The dust-to-gas ratio will rapidly 
reach a constant value, i.e., Zd — » yd as M s /M„ — » 00 if no dust 
destruction by SNe is present, and Z d — > yjy if it is. With the 
'secondary' dust included, Zd again approaches a constant value, 
this time given by 



Ji L 1 

V (JJV'z 1 

y \ y 



yd + eyz, 
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0, — -> 0. 
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To sum up, Zd has (within the present framework) an upper 
limit with only two exceptions: 

(i) a closed box without any dust destruction, 

(ii) a closed box with 'seconday' dust production. 

Since a real galaxy cannot strictly be a 'closed box', it is fair to 
assume that there is an upper limit to the gas-to-dust ratio Zd in 
most cases. 



2.6 Evolution of the dust-to-metals ratio 

A constant dust-to-metals ratio £ = Z A jZ, is commonly adopted in 
mode ls of galactic evolution (e.g. Edmunds & Ealesll998l ; |Pei et al.l 
1999). From observations, the dust-to-meta ls ratio £ is know n to be 
essentially constant in the local Universe dlssa et al] [T990). while 
a lower ( but simi l ar) ra ti o is found at somewhat h igher redshift 
(see, e.g. Iviadilol Il998l) . lEdmunds & Ealesl dl998h considered a 



simplistic mode l where f was fixed, but it was later shown by 
Edm undsl ( [200lh that £ could, at least in principle, vary significantly 
over time. However, the ^-evolution approaches a constant value in 
many cases, so for evolved galaxies of similar type one may expect 
an approximately universal £ . 

From Eqs. d 1 8125 1 > it is evident that the dust-to-gas ratio Z d for 
the different models can be expressed as functions of the metallicity 
Z. Hence, for a closed-box 



<T(Z) = 



1 - exp -v — 

yz 



zr 1 

v — + a>, 

yz 



yz 2yz' 



<5ism * 0, 



<5ism = 0, 



(27) 



and for the extreme infall case, 



?CZ) = 



yA 

- — he 
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Jz 



yz 
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+ w, Sism * 0, 



= 0. 



(28) 



From the equations above it is obvious that the evolution of f leads 
to different 'end-states' depending on the assumptions made for 
each model, f may thus reach a constant value, but can also grow 
without bound or approach zero as a galaxy evolves. 

Since very high dust-to-metals ratio s are not obser ved at any 
redshift (but see, e.g. lGalliano et al.ll2003Hviadilol20ol for exam- 
ples of variations), it seems the closed-box case with (5 1SM = and 
'secondary' non-stellar dust production is less likely, or star for- 
mation has to essentially cease at some point. It also appears that 
either the dust yield must be a significan fraction of the metal yield, 
or non-stellar 'secondary' dust production must play an important 
role in the local as well as the high-z Universe. For the Milky Way 
£mw = 0.3 - 0.5, which is roughl y the ratio foun d in most local 
galaxies of similar type (see, e.g. Ilssa et al1ll990h, and DLAs at 
z = 0.7 - 2.8 also show similar ratios (Vladilo 1998). 



3 RESULTS AND DISCUSSION 

Valiante et al. have suggested that stellar production of dust 

(in SNe and AGB stars) is high enough to explain the dust-mass 
estimates of quasar host galaxies at high z, in particular the host 
galaxy of SPSS Jl 148+525. But is this true under all circum- 
stances? Michalowski et alj ( feOlObl) have shown that the minimal 
yield needed to explain the observationally inferred dust masses in 
sub-millimeter galaxies at z > 4 is 15 - 65M Q per SN if the ob- 
servational constraints on the stellar and gas mass are assumed to 
be correct. Their result is marginally consistent wi th the theoret- 
ical maximum obtained by iTodini & Ferraral j20oTh . but does not 
account for any type of dust destruction. 

The main result of the present paper is that only the most ex- 
treme, but also less likely, scenarios for cosmic dust production 
can reproduce the dust masses inferred from observations of high-z 
galaxies. That is a scenario which requires essentially no dust de- 
struction - neither due to the reverse shock in SN-explosions nor 
in the ISM due to kinetic energy injection by SNe - and extensive 
non-stellar dust production seems necessary. 

Another important result is that galaxy formation by merger 
events (and accretion of pristine gas), i.e., a cosmologically mo- 
tivated rate of infall, slows down the build-up of the interstellar 
dust component quite significantly. This is partly due to the trivial 
fact that the ISM is diluted from the unenriched gas falling into the 
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Figure 2. Left: Comparison between evol ved infall models with 'secon dary' dust production (but without dust destruction) and observationally inferred dust 
masses of high-z galaxies (data taken from Michalowski et al. 2010a b). The light-blue region (appears light-gray in printed version) is only accessable using 
the 'observed' yield. The darker blue/gray region can only be reached using the theoretical 'maximum' yield and the hatched area in the middle corresponds 
to the region accessable using the 'high' yield. The plotted gas masses are twice the molecular hydrogen gas mass to account for neutral gas. For comparison 
the estimated upper limits of the dust and gas masses of the host/surroundings of GRB 090423 ( z = 8.2,| Stanway et alj|201ll) are overplotted, as well as the 
derived dust and gas masses for the less distant sub-millimeter galaxy SMM J02399-0136 (z = 2.8 lGenzel et alj2003l). Right: Same models as in the left panel, 
but compared to local galaxies, with and without the inclusion of neutral hydrogen in the gas mass (data taken from Devereux & Young 1990). 



galaxy, but also the fact that an infall scenario typically leads to 
significantly shorter dust destruction time scales (see Fig. [5}. 



3.1 Comparison with observed properties of high-z galaxies 

It is reasonable to restrict a general comparison with observations 
to constraints derived for the infall scenario, since in the current 
paradaign for galaxy formation galaxies are formed through mass 
assembly. Furthermore, it can be assumed that dust rich high-z 
galaxies are, despite their redshifts, in a relatively evolved state 
where a large fraction of the baryonic mass is stars, i.e., the time 
scale for their assembly short. A good starting point is therefore 
Eq. l|26} in Sect. fUl\ 

Using Eq. i26l with <5jsm = one can obtain upper lim- 
its to the dust masses possible in an evolved system for a given 
dust yield and efficiency e of 'secondary' dust production. With 
< e < 1 - Vd/yz as the allowed range, Eq. ((26} constrains the 
possible regions in the dust/gas-plane for given yields. In Fig. [2] 
evolved infall models with 'secondary' dust production (but with- 
out dust destruction) are co mpared to observationally inferred dust 
masses o f high-z galaxies dMichalowski et al. 2010a b) and local 
galaxies jDevereux & Youndl 19900 . 

The light-blue region (grey in printed version) is only accessi- 
ble using the 'observed' yield and the darker blue/grey region can 
o nly be reached using the theoretical 'maximum' yield according 
to lTodini & Ferraral ( l200l|) . The overlapping (hatched) region is ac- 
cessible using either one of the 'high' and 'maximum' yield de- 
pending on the amount of 'secondary' dust production. The dust 
content of the galaxies with the highest dust-to-gas ratios can only 
be explained with the theoretical 'maximum' yield in combination 



with a relatively large e. This will be discussed in more detail be- 
low. Note that the gas mass for the high-z galaxies are taken to be 
twice the observed IL-mass to account for a neutral gas component, 
which cannot be detected in these obj ects (the fraction of neutral 
gas has been estimated to be < 0.5, see lMcGreer et al.l201ll) . 

As can be seen in Fig. [2] the dust-to-gas ratios of high-z galax- 
ies are quite remarkable. Half of the sample considered here have 
Z d > 0.02, i.e., larger than the total metallicity of the solar neigh- 
bourhood. One may argue that this fact implies gas masses which 
are about one order of magnitude larger than those obtained through 
observational est imates of the molecular gas (see also models by 
iGall et alj|2010bl) . However, the uncertainties associated with the 
gas mass are certainly not on an order of magnitude scale. In fact, 
the total gas mass is not likely to be more than roughly a factor 2 
larger if atomic hydrogen and helium is added, and the uncertainty 
in the conversion factor between CO and H 2 (a) can at most be a 
few times larger than the value commonly used for high-z objects. 
A value around a ~ 1 i s normally used for high-z quasars (see, 
e.g. lDownes & Solomon 199jl), whi le for local spiral galaxies a, is 
typically taken to be ~ 4.6 (see, e.g. lSolomon & Barrett! 199~ll) . For 
local starburst galaxies a ~ 1, however. 

It seems non-stellar dust production/growth may be required 
to reproduce the high dust-to-gas ratios. Dust may grow in dense 
molecular clouds in the ISM through accretion onto see d parti- 
cles ( most likely originating from stars) given enough time (Draine 
1990). This component may be significant, possibly a more impor- 
tant contri bution than dust grown in stellar a tmospheres and SN 
remnamts dZhukovska. Gail & Trieloffl 120081 ; lAsano et al.l l201fj|) . 
Furthermore, lElvis et al. I d2002l) have suggested dust may be nu- 
cleated (and grow) in broad emission line clouds in outflowing 
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Figure 3. Star-formation history in the host galaxy of SDSS Jl 148+5251 
as predicted by the simulation of Li et al. (2007). The verical, dashed line 
marks z = 6.42 and the horizontal, dotted line shows the star-formation 
rate predicted by the Kennicutt-Schmidt law adopting a total gas mass of 
~ 10 10 MO (Walter et al. 2004). 



quasar winds. Both processes (growth in the ISM and formation 
in quasar winds) can be indirectly associated with star formation. 
Stars form primarely in dense molecular clouds which is also where 
dust grains in the ISM may grow through accretion. Quasar-wind 
outflow rates can be considered roug hly proportional t o the mass 
of the supermassive black hole (see lPipino et aljfeoi ll references 
therein) and are thus linked to the star formation rate (or growth 
of the stellar mass) assumi ng that the black hole is growing due t o 
the galactic mass assembly (DeBuhr et alj2010llSarria et al ■l20ld) . 
Hence, the non-stellar dust production/growth may effectively be 
modelled by introducing a modified (metallicity-dependent) yield 
as in Eq. H9\ . The possibility of grains forming in quasar ourflows 
and grains growing in the ISM should be investigated further by 
more detailed modelling before one concludes that this is the solu- 
tion to the dust formation problem in the early Universe, however. 

Although significant non-stellar dust production is an attrac- 
tive idea, much of the dust production problem at high z can be 
resolved by assuming the true gas mass may be a few times higher 
due to a significant atomic component and uncertainties in the 
derivation of the molecular component (see above), combined with 
the fact that dust masses are hard to quantify with high precision. 
It may very well be that the dust masses derived for, e.g., high-z 
quasars are too large due to the fact that the model for the dust 
emission is quite rudimentary. Dust mass estima tes are generally 
uncertain up to a factor of a few dSilva et al.ll998l) . Hence, the true 
dust-to-gas ratio Z d could, at least in principle, be an order of mag- 
nitude lower, which could then be explained by stellar dust produc- 
tion alone. If Z d is only a few times 1CT 3 , the required dust yield is 
reasonable from a theoretical point of view (c.f. the effect of reverse 
shocks) and consistent with results for local galaxies (see Fig.|2](. 
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Table 3. Baryonic mass components and IMFs used to model SDSS 
Jl 148+525. 
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4.3 
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0, 100, 200, 1000 
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91-120 
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1.6 
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0, 100, 200, 1000 


Normal 


121-150 


2.9 


■ 10 10 


1.6 


10 10 


0, 20, 100, 200 
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3.2 SDSS J1148+525 

The host galaxy of the SDSS Jl 148+5251 quasar has received 
considerable attention in the literature since it appears to have an 
unexpectedly large dust mass for its high redshift (z = 6.42), 
which seems to require a rather dramatic early star formation 
history if dust is mainly produced in stars (see, e.g., the mod- 
els by lDweketalJl2007l : IValiante et all 120091 ; iGall et alj l2010bl : 
|Pipinoetal.ll201ll) . This quasar host galaxy is used as a test bench 
also in the present paper. SDSS Jl 148+525 provides insight on the 
evolutionary states of dust-rich quasar hosts at high redshifts. 

Massive high-z galaxies seem to have reached a relatively 
evolvled state, despite their high redshifts, indicating that the y form 
in a ra pid series of mergers/infall events as in the model bv lLi et all 
(2007). The time scale of the build-up of the stellar mass is clearly 
not more than a few Myr. As a consequence, they have typically 
reached solar, or even supersolar, metallicities and the subsequent 
evolution (from the observed state until present time) will neither 
change the metallicity, nor the dust-to-gas ratio dramatically. 

The age of the Universe at z = 6.42 in the standard ACDM 
model is f a « 840 Myr and galaxy formation is here assumed to 
start at z = 14 corresponding to f a as 290 Myr. The host galaxy 
of SDSS J 1148+525 1 has in such c ase evolved for 550 Myr prior 
to its observed state. Ili et alj d2007l) constructed a merger tree, ex- 
tracted from a cosmological simulation, to follow the hierarchical 
mass assembly of a galaxy thought to be similar to t he host galaxy 
of SP SS Jl 148+5251. Here, just as in the work by IValiante et all 
(2009), their resultant star-formation history (SFH) is used as input 
(see Fig. |3). In practice, however, the quantity used is the time- 
integrated star-formation rate, i.e., the stellar mass, 

M s (t) = a f M t {t')dt', t = t a (z) - f a (z = 14), (29) 
Jo 

which was computed numerically from the d ata shown in Fi g.[3]us- 
ing a simple Romberg routine. In contrast to lValiante etal. (2009), 
a scaled-down version of this SFH is considered to meet observa- 
tional constraints on the total mass. 

Using the analytic model described in Sect. 12. 3| 150 param- 
eter configurations are considered (see Appendix |A) in order to 
find the most likely models given some ob servational constra ints to 
the dust mass in SDSS Jl 148+5251 (see I Valiante et alj|2009l. and 
references therein) and the dust-to-metals rati os found in da mped 
Lyo- absorbers at lower redshifts (1 < z < 2.5 Vladilo 1998J). Al- 
though inconsistent with a merger scenario, t he closed-box case is 
here considered b ecause (1) previous studies (IValiante et alj|2009t 
IGall et al.|[20T0bh have considered the SDSS Jl 148+5251 quasar 
host galaxy to be a 'closed box', (2) it is important to analyse the 
net differeces between an infall and closed-box scenario, since the 
evolution of the gas mass affects net dust-destruction efficiency (see 
Eq. [5} and (3) infall always lowers the efficiency of metal produc- 
tion TEdmundsIl 1990h . 

Only a few of the 150 models seem to reproduce both the 
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Figure 4. Best-fit models for SDSS Jl 148+525. The left panel shows the evolution of the dust mass. The grey region indicate the range of dust masses 
considered to be consistent with observations. The red circle with error bars show the average observationally deduced dust mass. The right panel shows 
the dust-to-metals ratio for the same models. The grey region again indicate the expected range of dust-to-metals ratios based on observations of DLAs at 
moderately high redshift. 



expected dust mass and a reasonable dust-to-metals ratio f (the 
models presented in Fig. |4). A plausibl e dust-mass range f or the 
SDSS Jl 148+525 host is that adopted bv lValiante etai] ( l2009l) . i.e., 
2-1 ■ 10 8 M o . Given the spread an d uncertainty of the d ust masses 
obtained by different authors (see IValiante et al.1l2009L and refer- 
ences therein), a dust mass in the range 10 8 ~ 9 M o is here considered 
to be consistent with the observations. A 'reasonable' range for the 
dust-to-metals ratio is here taken to be the sa me as for the s am- 
ple of DLAs at medium high redshift studied by Vladik] dl998l) . as 
£ is not changing much after the intial phase of galaxy evolution 
unles s the grain-growth rate is particularly low (see, e.g. lEdmundsl 
l200lL Fig. 4). Assumi ng dust of Galactic type, the DLAs consid- 
ered bv lVladild Jl998h places this quantity between 42% and 89% 
(typically 60%) of the Galactic v alue, where the latter is C, ~ 0.5 
dMever et al.ll 1 9981 : 1 Whitte^l 1 99 lh . 

The grid of models can be divided into five sub-sets of models 
based on the assumptions made regarding total mass, stellar mass, 
gas mass (at z = 6.42), dust destruction and the IMF (see Tables 
[T]and[3}. The first sub-set (models 1-30) represents an attempt to 
reprouce the metals and dust assuming the star formation b uilds up 
to rou ghly the same stellar mass as in the simulation by iLi et al.l 
d2007l) , while still maintaining the lower-limit gas mass suggested 
by observations. With model 1 1 as the only exception, non-stellar 
dust production appears necessary to reach above 10 8 M o of dust at 
z = 6.42. A few models (6, 12 and 26) do reach above M A = 1O 8 M 
and may be considered marginally consistent with observations in 
that regard. Model 26 is the only infall model without non-stellar 
dust production (among models 1-30) which is consistent with the 
dust mass derived from observations (see Fig. |4}. It is also consis- 
tent with the range of dust-to-metals ratios one would expect for a 
galaxy like SDSS Jl 148+525 (see Fig.|^ 

Changing the IMF from a normal iLarsonl i ll 9981) IMF (m c = 
0.3M o ) to at top-heavy ditto (ra c = 10M o ) increases both the 



amount of dust production as well as the dust destruction by SN 
shocks, as can be seen from models 31-60. Models 36, 42, 51 and 
57 are all reproducing the expected range for the dust mass, while 
all models with higher dust destruction efficiency (misM > 20M o ) 
predict too much dust destruction due to the significantly increased 
SN-rate. The resultant dust masses are in those cases at most a 
few times 10 7 M o . Models 31-60 are essentially ruled out because 
the resultant metallicities are too high to be realistic (Z = 1.47 
at z = 6.42) and inconsistent with the expected range of dust-to- 
metals ratios (see Fig. [6}. 

Models 61-90 ( see Fig. [6] and Table |A3ll are analogous to the 
numerical model by Valian te et al.l d2009l) in that they have similar 
stellar and gas masses (a total baryon mass of M tot = 1.3 • 10 12 M o ) 
and represent a closed-box scenario. The stellar dust yield re- 
ferred to as the 'hi gh' yield (see Table [2} and corresponds to the 
model preffered by IValiante et al.l d2009h . The resultant gas mass 
(M gas - 4.3 • 10"M o ) is sig nificantly higher than the mass inferred 
from observations (see, e.g. Walter et al. 2 004b. With this high g as 
mass, even the low ('observed') dust yield (see lGall et al . 2010b) is 
sufficient to reproduce the observed dust mass, and the 'maximal' 
yield leads to an over-production of dust in several cases (mod- 
els 71-73 and 86-88 ). This is the main reason why the model by 
IValiante et ail (2009) is marginally consistent with the high dust 
mass derived from observations without invoking dust growth in 
the ISM or any other seconday dust source. More gas simply means 
more dust without resulting in unrealistic dust-to-gas ratios. But as 
one can see in Fig. [6] the total amount of metals is generally too 
high to have dust-to-metals ratios in the expected range. 

Applying observational constraints on both the gas mass and 
the stellar mass (models 91-120) results in two plausible models. 
Both are models without dust destruction in the ISM and the 'max- 
imal' stellar yield (see Table |2j. Adding non-stellar dust produc- 
tion (as it is implemeted here) means no actual improvement for 
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the models with lower yields. It essentially requires that all metals 
in the ISM are turned into dust. 

Adding a top-heavy IMF (models 121-150, see Fig, IA5t leads 
to a result very similar to that of models 31-60 (also with a top- 
heavy IMF), i.e., the over-all metal production is not sufficient to 
reproduce the expected range of dust to metals. 

3.3 Infall vs. closed box 

Among the seven models of SDSS Jl 148+525 that reproduce the 
expected dust mass and dust-to-metals ratio (see previous section), 
only two are infall models. One may note as we ll that all the 
close d-box models presented here (and th e model bvlValiante et al.l 
2009) are conceptually inconsistent, since iLi et al.l d20071) consider 
a galaxy being formed by a sequence of merger events, where the 
smaller progenitor galaxies are very metal poor and have had very 
little star formation (M g /M tot « 1). The infall models (16-30, 46- 
60, 76-90, 106-120 and 136-150) are th erefore more appro priate. 
This fact was indeed acknowledged by IValiante et alj d2009l) . al- 
though the exact consequences were not discussed. 

From the present study it is clear, however, that infall models 
are not preferred over closed-box models in that they do not re- 
produce the dust mass and/or overall metallicity better than closed- 
box models. Infact, a merger/infall scenario requires a higher stellar 
dust yield in order to produce the same dust mass, which is due to 
the effects of the dilution by pristine infall (see, e.g. Edmunds 1990) 
and the fact that the dust destruction time scale is proportional to 
gas mass. The latter is important during the early phase of evolu- 
tion, since in a closed-box scenario the initial gas mass is always 
equal to the total mass of baryons. The gas mass is never as high 
as during the initial phase of the gas mass evolution, which means 
that the dust destruction time scale is significantly longer at early 
times (when the rate of dust production is at its highest) compared 
to the infall scenario (see Eq.|5]and the example in Fig. [5}- 

In the model of SDSS Jl 148+525, the intense star formation 
expected during the first 400 Myr would produce a lot of stel- 
lar dust, while T d would still be relatively long on average in the 
closed-box case. As a consequence, the effective dust production 
have to be higher. If an infall scenario is adopted, the gas mass 
during the early evolution is lower even if the final total mass of 
baryons is the same. For example, in models 1-30 the dust destruc- 
tion time scale is about 15 times shorter for infall models compared 
to closed-box models and therefore infall models have a lower dust- 
formation efficiency. 

Infall models always require a larger metal yield than the 
closed-box models in order to contain the same amount metals. As 
discuss ed above, the d ust-to-metals ratio £ should be roughly that of 
DLAs dVladilol 19981) . The lower net-production of dust for a given 
stellar mass of infall models tends to result in higher values (see 
Fig. [6]l. The infall models presented here are perhaps unrealistic in 
that they assume all gas ejected by evolved stars and all infalling 
gas is converted into new stars (keeping dM g /dt = 0). But the dust 
yield required to produce the estimated dust mass is at the high 
end and the star-formation efficiency will under all circumstances 
be high in quasar host galaxies at high redshifts. 



3.4 Quasar outflows 

As argued in Sect. 12.21 the presence of a galactic wind driven by 
SNe will have roughly the same effect on the dust mass as dust de- 
struction by SNe has, thus making it even harder to obtain the high 
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Figure 5. Dust destruction time scale as function of time for a closed-box 
and an extreme infall scenario using the star-formation history resulting 
from the simulation of Li et al. (2007). 

dust masses in some high-z galaxies. However, the outfl ows associ- 
ated w ith quasars may also have a totally different effect. lElvis et al.l 
d2002h suggested that dust may be created in broad emission line 
clouds associated with quasar outflows. As opposed to the gas in 
the ISM, these gas clouds could in principle (remains to be proven) 
have both nucleation and growth of dust grains. This fact makes 
quasar outflows a possible source for the 'missing dust' in models 
of dust formation in high-z quasar host galaxies. The models pre- 
sented here suggest the efficiency of dust production must be quite 
high, i.e., a large fraction of the atomic metals in the gas must be 
converted into dust (essentially all the metals in the galaxy need to 
be in dust in many cases, see tables in Appendix lAl>, This might be 
a problem, since one can only have as much metals as can be pro- 
duced by the stars and the total gas mass in a quasar outflow is not 
likely to be very large. Thus quasar outflows may contribute, but 
the contri bution can be con sidered negligible (see also modelling 
results bv |Pipinoetai1l201ll) . 

3.5 A top-heavy IMF? 

There is a variety of evidence for top-heavy IMFs in starburst 
enviro nments (see, e.g.. |Pab ringhausen et al. 20091). IValiante et all 
d2009h showed that in their model LIM stars contributed about 80% 
of the dust, using a iLarsonl ( fl998h IMF, which is similiar to th at of 
ISalpetej [l955) except at low masses where the lLarsonl d 1998b IMF 
has a smooth turn-over. But with a top-heavy IMF the situation is 
quite different. 

From observations it seems LIM stars are producing signif- 
icantly more dust per stellar mass than HM stars (cf. the 'ob- 
served' yield in Table [2}- Hence, a top-heavy IMF would under 
these circumstances result in a lower total stellar dust yield y A . If 
the dust production is dominated by HM stars, on the other hand, 
a top-heavy IMF leads to a significantly increased dust production 
dGalletalfc OlOb). However, a larger fraction of HM stars also in- 
creases the SN-rate, which in turn leads to more interstellar dust 
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destruction due to SN-shocks (see Sect.JO] Fig.|6] Models 31-60 
and 121-150 in Appendix [At. But if this type of dust destruction is 
negligible, an initially top-heavy IMF may help to boost dust pro- 
duction in the early Universe. In general, though, a top-heavy IMF 
means no particular improvment since (1) dust destruction due SN- 
shocks is expected and (2) the dust-to-metals ratio changes in an 
unfavourable direction (again, see Fig. [6] where it is clear that the 
ratios are too low). 



4 SUMMARY AND CONCLUSIONS 

Stars can only produce a limited amount of dust and the time-scale 
of survival of dust in the ISM may not be long enough. Thus, there 
are three fundamental reasons why stellar sources probably cannot 
explain the dust masses observed at several high-z objects: 

(i) Stellar yields must under most circumstances be close to the 
theoretical maximum for simple dust evolution models to be con- 
sistent with dust masses derived from observations. 

(ii) A top-heavy IMF means no paricular improvement since 
the SN-rate increases and therefore the rate of dust destruction in- 
creases as well. Consequently the dust-to-metals ratio changes in 
an unfavourable fashion. 

(iii) Infall of gas (which is more likely than a closed-box sce- 
nario) may shorten the dust-destruction time scale and therby lower 
the effective stellar dust production. 

The simple analytical models presented in this paper suggest the 
dust-mass problem in high-z galaxies cannot be solved by only 
adding the co ntribution from AGB stars. The inclusion of AGB 
stars by made Valiante et al. (2009) adds new and important infor- 
mation, but it does not provide a complete and viable solution to 
the dust product i on pro blem at high redshifts. 

lElvis et all d2002l) have suggested dust may be created in 
broad emission line clouds in an outflowing quasar wind. One may 
also assume dust grains can form, or at least grow, in the interstel- 
lar gas of quasar host galaxies. However, neither hypotheses can 
properly solve the dust production problem, because the amount of 
metals needed as raw-material for dust formation cannot be pro- 
duced by stars alone in many cases. 

There is essentially only two ways out of the dilemma de- 
scribed above, is either significant non-stellar dust production (e.g., 
growth of dust grains in the ISM or dust formation associated with 
quasar outflows) or signifiant systematic errors in the dust and 
(molecular) gas masses derived from observations. In a few cases 
it may be sufficient to assume no dust destruction, but this cannot 
solve the problem in general. 
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Table Al. Parameter values for models with a normal IMF and observed gas mass. Masses are given for z = 6.42. 



Model 


a 


e 


>'d 






yz 




"»ISM Wo\ 


<5ism 


V 


M g , 


[M G ] 


M tot 


[M Q ] 


Description 


1 


0.63 




8.81 


10" 


-4 


3.00 


1 a-2 

10 





0.00 




1.6 


10 


8.85 


10 


Closed box, no dust destruction. 


2 


0.63 




8.81 


10" 


-4 


3.00 


1 a-2 

10 


100 


1.00 


1.59 


1.6 


i aIO 
10 


8.85 


1 a1 l 
10 


Closed box. 


3 


0.63 




8.81 


10" 


-4 


3.00 


10 


200 


2.00 


3.17 


1.6 


i aIO 
10 


8.85 


1 a1 l 
10 


Closed box. 


4 


63 




8 81 


10" 


-4 


3 00 


IQ-2 


1 000 


10 


15 9 


1 6 


1()10 


o.o J 


in" 


Closed box. 


5 


0.63 


0.35 


8.81 


10" 


-4 


3.00 


10- 2 





0.00 




1.6 


10 10 


8.85 


10 u 


Closed box, 'secondary dust'. 


6 


0.63 




1.52 


10" 


-3 


3.00 


1 A— 2 

10 





0.00 




1.6 


10 


8.85 


1 All 

10 


Closed box, no dust destruction. 


7 


0.63 




1.52 


10" 


-3 


3.00 


1 a-2 

10 


100 


1.00 


1.59 


1.6 


i aIO 
10 


8.85 


1 All 
10 


Closed box. 


8 


0.63 




1.52 


10" 


-3 


3.00 


i a-2 
10 


200 


2.00 


3.17 


1.6 


10 


8.85 


1 AH 
10 


Closed box. 


9 


0.63 




1.52 


10" 


-3 


3.00 


10~ 2 


1000 


10.0 


15.9 


1.6 


10 10 


8.85 


10 u 




10 


0.63 


0.30 


1.52 


10" 


-3 


3.00 


io- 2 





0.00 




1.6 


10 10 


8.85 


10 u 


Closed box, 'secondary dust'. 


11 


0.63 




1.08 


10" 


-2 


3.00 


1 a-2 

10 





0.00 




1.6 


10 


8.85 


1 All 

10 


Closed box, no dust destruction. 


12 


0.63 




1.08 


10" 


-2 


3.00 


1 a-2 
10 


100 


1.00 


1.59 


1.6 


10 


8.85 


1 All 

10 


Closed box. 


13 


0.63 




1.08 


10" 


-2 


3.00 


1 a-2 
10 


200 


2.00 


3.17 


1.6 


10 


8.85 


1 All 

10 


Closed box. 


14 


0.63 




1.08 


10" 


-2 


3.00 


10~ 2 


1000 


10.0 


15.9 


1.6 


10 10 


8.85 


10" 


C~ 'IacpH hnv 


15 


0.63 


0.30 


1.08 


10" 


-2 


3.00 


io- 2 


100 


1.00 


1.59 


1.6 


10 10 


8.85 


10" 


Closed box, 'secondary dust'. 


16 


0.63 




8.81 


10" 


-4 


3.00 


1 a-2 

10 





0.00 




1.6 


1 aIO 

10 


8.85 


1 All 

10 


Mali, no dust destruction. 


17 


0.63 




8.81 


10" 


-4 


3.00 


1 a-2 

10 


100 


1.00 


1.59 


1.6 


10 


8.85 


1 All 

10 


Infall 


1 8 


0.63 




8.81 


10" 


-4 


3.00 


1 a-2 
10 


200 


2.00 


3.17 


1.6 


mlO 
10 


8.85 


1 a1 1 
10 


Infall 


J J 


A f^X 

U.O J 




O.O 1 


1 u 


-4 


1 AA 


i a-2 

1U 


1 AAA 


1 A A 
1U.U 


1 ^ Q 


1 fi 




o.o J 


i a1 1 
i\J 


Tnfiill 

iniaii 


20 


0.63 


0.85 


8.81 


10" 


-4 


3.00 


io- 2 





0.00 




1.6 


io 10 


8.85 


10" 


Infall, 'secondary dust' . 


21 


0.63 




1.52 


10" 


-3 


3.00 


io- 2 





0.00 




1.6 


io 10 


8.85 


10" 


Infall, no dust destruction. 


22 


0.63 




1.52 


io- 


-3 


3.00 


io- 2 


100 


1.00 


1.59 


1.6 


10 10 


8.85 


10" 


Infall 


23 


0.63 




1.52 


io- 


-3 


3.00 


io- 2 


200 


2.00 


3.17 


1.6 


10 10 


8.85 


10" 


Infall 


24 


0.63 




1.52 


10" 


■3 


3.00 


10~ 2 


1000 


10.0 


15.9 


1.6 


10 10 


8.85 


10" 


Infall 


25 


0.63 


0.80 


1.52 


10" 


-3 


3.00 


io- 2 





0.00 




1.6 


10 10 


8.85 


10" 


Infall, 'secondary dust'. 


26 


0.63 




1.08 


10" 


-2 


3.00 


io- 2 





0.00 




1.6 


10 10 


8.85 


10" 


Infall, no dust destruction. 


27 


0.63 




1.08 


10" 


-2 


3.00 


io- 2 


100 


1.00 


1.59 


1.6 


10 10 


8.85 


10" 


Infall 


28 


0.63 




1.08 


io- 


-2 


3.00 


io- 2 


200 


2.00 


3.17 


1.6 


10 10 


8.85 


10" 


Infall 


29 


0.63 




1.08 


10" 


-2 


3.00 


10~ 2 


1000 


10.0 


15.9 


1.6 


IO 10 


8.85 


10" 


Infall 


30 


0.63 


0.50 


1.08 


io- 


-2 


3.00 


io- 2 





0.00 




1.6 


IO 10 


8.85 


10" 


Infall, 'secondary dust'. 
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Figure Al. Models using the star-formation history resulting from the simulation of Li et al. (2007). The filled circle with error bars show the estimated dust 
mass of SDSSJ1 148+525. 
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Table A2. Parameter values for models with a top-heavy IMF and observed gas mass. Masses are given for z = 6.42. 



Model 


a 


e 


yd 






yz 




'"ISM [Mq] 


<5ism 


V 


Mgas [M Q ] 


M tot 


[M ] 


Description 


31 


0.10 




3.51 


10" 


-3 


6.50 


10 





0.00 




1.6 


10 10 


1.54 


mil 
10 


Closed box, no dust destruction. 


32 


0.10 




3.51 


10" 


-3 


6.50 


10 _1 


20 


1.32 


1.59 


1.6 


10 10 


1.54 


1 All 

10 


Closed box. 


33 


0.10 




3.51 


10" 


-3 


6.50 


10 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


inll 
10 


Closed box. 


34 


10 




3 51 


10" 


-3 


6 50 




200 


13 2 


15.9 


1.6 ■ 


10 


1.54 


10" 


v^ltJftCU UtJA. 


35 


0.10 


0.25 


3.51 


10" 


-3 


6.50 


10 1 


20 


1.32 




1.6 


i r\10 

10 


1.54 


10" 


Closed box, 'secondary dust'. 


36 


0.10 




1.55 


10" 


-2 


6.50 


10" 1 





0.00 




1.6 


10 10 


1.54 


i nil 
10 


Closed box, no dust destruction. 


37 


0.10 




1.55 


10" 


-2 


6.50 


10 


20 


1.32 


1.59 


1.6 


10 10 


1.54 


inll 
10 


Closed box. 


38 


0.10 




1.55 


10" 


-2 


6.50 


10 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


mil 
10 


Closed box. 


39 


10 




1 55 


10" 


-2 


6 50 


10 -1 


200 


13 2 


15.9 


1.6 


i r\10 

10 


1.54 


10" 


UIUlICU UtJA. 


40 


0.10 


0.25 


1.55 


10" 


-2 


6.50 


10~' 


20 


1.32 




1.6 


i nlO 
10 


1.54 


10" 


Closed box, 'secondary dust'. 


41 


0.10 




1.88 


10" 


-1 


6.50 


10 





0.00 




1.6 


10 10 


1.54 


1 nl 1 
10 


Closed box, no dust destruction. 




0.10 




1 QQ 

1.88 


10" 


-1 


6.50 


10" 




1.32 


1 £ o 

1.59 


1.6 


10 10 


1.54 


inll 
10 


Closed box. 


43 


0.10 




1.88 


10" 


-1 


6.50 


10 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


mil 
10 


Closed box. 


44 


0.10 




1.88 


10" 


-1 


6.50 


10" 1 


200 


13.2 


15.9 


1.6 


10 


1.54 


10" 




45 


0.10 


0.12 


1.88 


10" 


-1 


6.50 


10 ! 


20 


1.32 


1.59 


1.6 


10 


1.54 


10 n 


Closed box, 'secondary dust'. 


46 


0.10 




3.51 


10" 


-3 


6.50 


10 





0.00 




1.6 


10 10 


1.54 


i nil 
10 


Infall, no dust destruction. 


47 


0.10 




3.51 


10" 


-3 


6.50 


10 


20 


1.32 


1.59 


1.6 


10 10 


1.54 


mil 
10 


Infall 


48 


0.10 




3.51 


10" 


-3 


6.50 


10 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


mil 
10 


Infall 


4V 


n i n 




^ ^1 
j.j i 


1 IV 
LU 


-3 


A ^n 


1 n~l 


4XJO 


1 1 T 


1 ^ Q 


1.6 


10 10 


1 

I.J4 


inn 


inran 


50 


0.10 


0.55 


3.51 


10" 


-3 


6.50 


IQ- 1 


20 


1.32 




1.6 


10 10 


1.54 


10" 


Infall, 'secondary dust'. 


51 


0.10 




1.55 


10" 


-2 


6.50 


10 1 





0.00 




1.6 


10 10 


1.54 


10" 


Infall, no dust destruction. 


52 


0.10 




1.55 


10" 


-2 


6.50 


10~' 


20 


1.32 


1.59 


1.6 


10 10 


1.54 


10" 


Infall 


53 


0.10 




1.55 


10" 


-2 


6.50 


10~' 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


10" 


Infall 


54 


0.10 




1.55 


10" 


-2 


6.50 


10"' 


200 


13.2 


15.9 


1.6 


10 10 


1.54 


10" 


Infall 


55 


0.10 


0.55 


1.55 


10" 


-2 


6.50 


10"' 


20 


1.32 




1.6 


10 10 


1.54 


10" 


Infall, 'secondary dust'. 


56 


0.10 




1.88 


10" 


-1 


6.50 


10"' 





0.00 




1.6 


10 10 


1.54 


10" 


Infall, no dust destruction. 


57 


0.10 




1.88 


10" 


-1 


6.50 


10"' 


20 


1.32 


1.59 


1.6 


10 10 


1.54 


10" 


Infall 


58 


0.10 




1.88 


10" 


-1 


6.50 


10~' 


100 


6.60 


3.17 


1.6 


10 10 


1.54 


10" 


Infall 


59 


0.10 




1.88 


10" 


-1 


6.50 


10~' 


200 


13.2 


15.9 


1.6 


10 10 


1.54 


10" 


Infall 


60 


0.10 


0.30 


1.88 


10" 


-1 


6.50 


10"' 


20 


1.32 




1.6 


10 10 


1.54 


10" 


Infall, 'secondary dust'. 
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Figure A2. Same as Fig. lAll but for the models with a top-heavy IMF. 
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Table A3. Parameter values for models with a normal IMF and a total mass of 1.3 ■ 10 12 M o jLietalj2 007: Valian te et alj20O9h . Masses are given for z = 6.42. 



Model 


a 


e 


y<i 






yz 






rmsM [M Q ] 


<?ISM 


V 


M m [M e ] 


M lol 


[M ] 


Description 


61 


0.63 


- 


8.81 


10" 




3.00 


10" 







0.00 


- 


4.3 


10" 


1.30 


10 12 


Closed box, no dust destruction. 


62 


0.63 


- 


8.81 


10" 




3.00 


10" 




100 


1.00 


1.59 


4.3 


10" 


1.30 


10 12 


Closed box. 


63 


0.63 


- 


8.81 


10" 




3.00 


10" 




200 


2.00 


3.17 


4.3 


10" 


1.30 


10 12 


Closed box. 


64 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


1000 


10.0 


15.9 


4.3 


10" 


1.30 


10 


Closed box. 






05 


8 81 




-4 


3 on 


1 0* 


-2 


LOO 


1 00 


1 S9 


4.3 


10" 




m 12 


Closed box, secondciry dust . 


66 


0.63 


- 


1.52 


10" 




3.00 


10" 







0.00 




4.3 


10" 


1.30 


10 12 


Closed box, no dust destruction. 


67 


0.63 


- 


1.52 


10" 




3.00 


10" 




LOO 


1.00 


1.59 


4.3 


10" 


1.30 


10 12 


Closed box. 


68 


0.63 


■ 


1.52 


10" 


_3 


3.00 


10" 


_2 


200 


2.00 


3.17 


4.3 


10 11 


1.30 


10 12 


Closed box. 


/n 
69 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


1000 


10.0 




4.3 


10" 


1.30 


10 


Closed box. 


70 


63 


06 


1 52 


10" 


-3 


3 00 


10" 


-2 


200 


2 00 


3 17 


4.3 


10" 


1 30 


10 12 


L,1UliCU UtJA, ftCUtJllUdl V LILLAL . 


71 


0.63 


- 


1.08 


10" 




3.00 


10' 







0.00 


- 


4.3 


10" 


1.30 


10 12 


Closed box, no dust destruction. 


72 


0.63 


- 


1.08 


10" 




3.00 


10" 




100 


1.00 


1.59 


4.3 


10" 


1.30 


10 12 


Closed box. 


73 


0.63 




1.08 


10" 




3.00 


10" 




200 


2.00 


3.17 


4.3 


10" 


1.30 


10 12 


Closed box. 




0.63 




i no 
1 .06 


10" 


-2 


3.00 


1 0* 


-2 


1000 


10.0 


1 £ o 

15.9 


4.3 


10" 


1.30 


to 


Closed box. 


7S 




12 


1 08 




-2 


3 on 


1 0" 


-2 


1000 


10 


1 S 9 


4.3 


10" 


1 30 


10 12 


Closed box, secondary dust . 


76 


0.63 


- 


8.81 


10" 




3.00 


10" 







0.00 


- 


4.3 


10" 


1.30 


10 12 


Infall, no dust destruction. 


77 


0.63 


- 


8.81 


10" 




3.00 


10" 




100 


1.00 


1.59 


4.3 


10 11 


1.30 


10 12 


Infall 


78 


0.63 




8.81 


lo- 




3.00 


10" 




200 


2.00 


3.17 


4.3 • 


mil 


1.30 


10 12 


Infall 


79 


0.63 




8.81 


ur 


-4 


3.00 


10* 


-2 


1000 


10.0 


15.9 


4.3 


10" 


1.30 


10 12 


Infall 


80 


0.63 


0.12 


8.81 


10" 


-4 


3.00 


io- 


-2 


200 


2.00 


3.17 


4.3 


10" 


1.30 


10 12 


Infall, 'secondary dust'. 


81 


0.63 




1.52 


10" 


-3 


3.00 


ur 


-2 





0.00 




4.3 


10" 


1.30 


10 12 


Infall, no dust destruction. 


82 


0.63 




1.52 


10" 


-3 


3.00 


ur 


-2 


100 


1.00 


1.59 


4.3 


10" 


1.30 


10 12 


Infall 


83 


0.63 




1.52 


lo- 


-3 


3.00 


io- 


-2 


200 


2.00 


3.17 


4.3 


10 11 


1.30 


10 12 


Infall 


84 


0.63 




1.52 


ur 


-3 


3.00 


10" 


-2 


1000 


10.0 


15.9 


4.3 


10" 


1.30 


10 12 


Infall 


85 


0.63 


0.10 


1.52 


10" 


-3 


3.00 


10" 


-2 


200 


2.00 


3.17 


4.3 


10" 


1.30 


10 12 


Infall, 'secondary dust'. 


86 


0.63 




1.08 


ur 


-2 


3.00 


io- 


-2 





0.00 




4.3 


10" 


1.30 


10 12 


Infall, no dust destruction. 


87 


0.63 




1.08 


io- 


-2 


3.00 


io- 


-2 


100 


1.00 


1.59 


4.3 


10" 


1.30 


10 12 


Infall 


88 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


200 


2.00 


3.17 


4.3 


10" 


1.30 


10 12 


Infall 


89 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


1000 


10.0 


15.9 


4.3 


10" 


1.30 


10 12 


Infall 


90 


0.63 


0.20 


1.08 


10" 


-2 


3.00 


10" 


-2 


1000 


10.0 


15.9 


4.3 


10" 


1.30 


10 12 


Infall, 'secondary dust'. 
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Figure A3. Same as Fig.ET] but for the models with a normal IMF and a total mass of 1.3 ■ 10 12 M G ihi et alj2 007; Valian te et alj2009h . 
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Table A4. Parameter values for models with a normal IMF and dynamically constrained total mass. Masses are given for z = 6.42. 



Model 


a 


€ 


yd 






yz 






"1ISM [Mq] 


<5lSM 


V 




[M G ] 


M lot 


[M Q ] 


Description 


9 1 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 





0.00 




1.6 ■ 


10 


4.5 • 


inlO 
10 


Closed box, no dust destruction. 


92 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


100 


1.00 


1.59 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 

10 


Closed box. 


93 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


94 


0.63 




8 81 


10" 


-4 


3 00 


10" 


-2 


1000 


10 


15 9 


1.6 ■ 


10 10 


4.5 ■ 


10 10 


V^lUoCU UtJA. 


95 


0.63 


1.50 


8.81 


10" 


-4 


3.00 


10" 


-2 





0.00 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


96 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 





0.00 




1.6 • 


i nlO 
10 


4.5 ■ 


inlO 
10 


Closed box, no dust destruction. 


97 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


100 


1.00 


1.59 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


98 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 ■ 


10 


4.5 ■ 


inlO 
10 


Closed box. 


99 


63 




1.52 


10" 


-3 


3 00 


10" 


-2 


1000 


10 


15 9 


1.6 ■ 


10 10 


4 5 ■ 


10 10 


v^lvJJiCU UUA. 


100 


0.63 


1.50 


1.52 


10" 


-3 


3.00 


10" 


-2 


() 


0.00 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


101 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 





0.00 




1.6 ■ 


inlO 
10 


4.5 • 


inlO 
10 


Closed box, no dust destruction. 




0.63 




1 .08 


10" 


-2 


3.00 


10" 


-2 


100 


1.00 


1.59 


1 

1.6 ■ 


inlO 
1U 


4.5 - 


inlO 
1U 


Closed box. 


103 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


104 


0.63 




1 08 


10" 


-2 


3 00 


10" 


-2 


1000 


10 


15 9 


1.6 ■ 


10 10 


4.5 • 


10 10 


i^iuacu uoa. 


105 


0.63 


1.00 


1.08 


10" 


-2 


3.00 


10" 


-2 





0.00 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


106 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 





0.00 




1.6 ■ 


inlO 
10 


4.5 • 


inlO 
10 


Infall, no dust destruction. 


107 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


100 


1.00 


1.59 


1.6 ■ 


inlO 
10 


4.5 • 


inlO 
10 


Infall 


108 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Infall 


109 


0.63 




8.81 


10" 


-4 


3.00 


10" 


-2 


1000 


10.0 


15.9 


1.6- 


10 10 


4.5- 


10 10 


Infall 


110 


0.63 


0.97 


8.81 


10" 


-4 


3.00 


10" 


-2 





0.00 




1.6 • 


10 10 


4.5- 


10 10 


Infall, 'sec. dust', max. e. 


111 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 





0.00 




1.6- 


10 10 


4.5 ■ 


10 10 


Infall, no dust destruction. 


112 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


100 


1.00 


1.59 


1.6- 


10 10 


4.5- 


10 10 


Infall 


113 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 ■ 


10 10 


4.5- 


10 10 


Infall 


114 


0.63 




1.52 


10" 


-3 


3.00 


10" 


-2 


1000 


10.0 


15.9 


1.6- 


10 10 


4.5- 


10 10 


Infall 


115 


0.63 


0.95 


1.52 


10" 


-3 


3.00 


10" 


-2 





0.00 




1.6 • 


10 10 


4.5- 


10 10 


Infall, 'sec. dust', max. e. 


116 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 





0.00 




1.6 ■ 


10 10 


4.5- 


10 10 


Infall, no dust destruction. 


117 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


100 


1.00 


1.59 


1.6 • 


10 10 


4.5- 


10 10 


Infall 


118 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


200 


2.00 


3.17 


1.6 • 


10 10 


4.5- 


10 10 


Infall 


119 


0.63 




1.08 


10" 


-2 


3.00 


10" 


-2 


1000 


10.0 


15.9 


1.6- 


10 10 


4.5- 


10 10 


Infall 


120 


0.63 


0.64 


1.08 


10" 


-2 


3.00 


10" 


-2 





0.00 




1.6- 


10 10 


4.5- 


10 10 


Infall, 'sec. dust', max. 6. 
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Figure A4. Same as Fig. lAll but for the models with a normal IMF and a total mass of 4.5 ■ 10 10 M Q as obtained from dynamical estimates ^Walter etal .120041) . 
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Table A5. Parameter values for models with a top-heavy IMF and dynamically constrained total mass. Masses are given for z = 6.42. 



Model 


a 


e 








yz 




"1ISM [M Q ] 


i5ism 


V 


M ga 


[M G ] 


M to , 


[M Q ] 


Description 


121 


0.10 




8.81 


10" 


-4 


6.50 


10 





0.00 




1.6 ■ 


inlO 
10 


4.5 • 


inlO 
10 


Closed box, no dust destruction. 


122 


0.10 




8.81 


10" 


-4 


6.50 


10 


20 


1.32 


1.59 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


123 


0.10 




8.81 


10" 


-4 


6.50 


10 


100 


6.60 


3.17 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


124 


10 




8 81 


1 0" 


-4 


6 50 




9fYl 
— wu 




1 S Q 
l J.y 


1 6 ■ 


10 10 




10 10 


Closed box. 


125 


0.10 


0.55 


8.81 


10" 


-4 


6.50 


io~' 


20 


1.32 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


126 


0.10 




1.52 


10" 


-3 


6.50 


10 





0.00 




1.6 ■ 


10 


4.5 • 


inlO 
10 


Closed box, no dust destruction. 


127 


0.10 




1.52 


10" 


-3 


6.50 


10 


20 


1.32 


1.59 


1.6 ■ 


inlO 
10 


4.5 ■ 


10 


Closed box. 


128 


0.10 




1.52 


10" 


-3 


6.50 


10 


100 


6.60 


3.17 


1 .6 ■ 


inlO 
10 


4.5 • 


10 


Closed box. 


129 


10 




1.52 


10" 


-3 


6 50 




200 


13 2 


15.9 


1.6 ■ 


10 10 


4 5 ■ 


10 10 


v^lvJJiCU UUA. 


130 


0.10 


0.50 


1.52 


10" 


-3 


6.50 


ur 1 


20 


1.32 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


131 


0.10 




1.08 


10" 


-2 


6.50 


10 





0.00 




1.6 ■ 


mlO 
10 


4.5 ■ 


inlO 
10 


Closed box, no dust destruction. 


132 


0. 10 




1 .08 


10" 


-2 


6.50 


10" 


20 


1.32 


1.59 


1.6 ■ 


1U 


4.5 - 


inlO 

tu 


Closed box. 


133 


0.10 




1.08 


10" 


-2 


6.50 


10 


100 


6.60 


3.17 


1.6 ■ 


inlO 
10 


4.5 ■ 


inlO 
10 


Closed box. 


134 


0.10 




1.08 


10" 


-2 


6.50 


10 _1 


200 


13.2 


15.9 


1.6 • 


10 10 


4.5 ■ 


10 10 




135 


0.10 


0.25 


1.08 


10" 


-2 


6.50 


10~' 


20 


1.32 


1.59 


1.6 ■ 


10 10 


4.5 ■ 


10 10 


Closed box, 'secondary dust'. 


136 


0.10 




8.81 


10" 


-4 


6.50 


10 _1 





0.00 




1.6 ■ 


mlO 
10 


4.5 ■ 


inlO 
10 


Infall, no dust destruction. 


137 


0.10 




8.81 


10" 


-4 


6.50 


10" 1 


20 


1.32 


1.59 


1.6 • 


10 


4.5 ■ 


inlO 
10 


Infall 


138 


0.10 




8.81 


10" 


-4 


6.50 


10 


100 


6.60 


3.17 


1.6 ■ 


mlO 
10 


4.5 ■ 


mlO 
10 


Infall 


139 


0.10 




8.81 


10" 


-4 


6.50 


10-' 


200 


13.2 


15.9 


1.6- 


10 10 


4.5- 


10 10 


Infall 


140 


0.10 


0.70 


8.81 


10" 


-4 


6.50 


ur 1 


20 


1.32 




1.6 • 


10 10 


4.5- 


10 10 


Infall, 'secondary dust' . 


141 


0.10 




1.52 


lo- 


-3 


6.50 


10-' 





0.00 




1.6 ■ 


10 10 


4.5 ■ 


10 10 


Infall, no dust destruction. 


142 


0.10 




1.52 


io- 


-3 


6.50 


10-' 


20 


1.32 


1.59 


1.6- 


10 10 


4.5 ■ 


10 10 


Infall 


143 


0.10 




1.52 


ur 


-3 


6.50 


10~' 


100 


6.60 


3.17 


1.6 ■ 


10 10 


4.5- 


10 10 


Infall 


144 


0.10 




1.52 


10" 


-3 


6.50 


io-' 


200 


13.2 


15.9 


1.6- 


10 10 


4.5- 


10 10 


Infall 


145 


0.10 


0.65 


1.52 


10" 


-3 


6.50 


lo- 1 


20 


1.32 




1.6 • 


10 10 


4.5- 


10 10 


Infall, 'secondary dust' . 


146 


0.10 




1.08 


10" 


-2 


6.50 


10-' 





0.00 




1.6 ■ 


10 10 


4.5- 


10 10 


Infall, no dust destruction. 


147 


0.10 




1.08 


10" 


-2 


6.50 


10"' 


20 


1.32 


1.59 


1.6 • 


10 10 


4.5 ■ 


10 10 


Infall 


148 


0.10 




1.08 


10" 


-2 


6.50 


ur 1 


100 


6.60 


3.17 


1.6 • 


10 10 


4.5 ■ 


10 10 


Infall 


149 


0.10 




1.08 


10" 


-2 


6.50 


lo- 1 


200 


13.2 


15.9 


1.6 ■ 


10 10 


4.5 ■ 


10 10 


Infall 


150 


0.10 


0.30 


1.08 


10" 


-2 


6.50 


lo- 1 


20 


1.32 




1.6- 


10 10 


4.5- 


10 10 


Infall, 'secondary dust' . 
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Figure A5. Same as Fig. lAll but for the models with a top4ieavy IMF and a total mass of 4.5 • 10 I0 M Q as obtained from dynamical estimates (Walter et al. 

|2004|) . 
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